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In creating his gravitational field equations Einstein assumed without justification that inertial 
mass, even in its equivalent form as energy, is a source of gravity. Denying that assumption allows 
modifying the field equations to a form in which a positive cosmological constant appears as a 
uniform density of gravitationally repulsive matter. Field equations with both positive and negative 
active gravitational mass densities of both primordial and continuously created matter incorporated 
along with two scalar fields coupled to geometry yield cosmological solutions that exhibit infiation, 
deceleration, coasting, acceleration, and a 'big bounce' instead of a 'big bang'. The repulsive matter 
is identified as the back sides of the 'drainholes' introduced by the author in 1973 as solutions of 
those same field equations. Drainholes are topological tunnels in space which gravitationally attract 
on their front, entrance sides and repel more strongly on their back, exit sides. The front sides serve 
in part as the unseen particles of 'dark matter' needed to hold together the large scale structures seen 
in the universe, the back sides, as the misnamed 'dark energy' driving the accelerating expansion 
of the universe. Formation of cosmic voids, walls, filaments, and nodes is attributed to expulsion of 
drainhole entrances from voids made and populated by drainhole exits, and subsequent aggregation 
of the entrances on the void boundaries. All of these cosmological entities are found wrapped in 
one coherent package, namely, the field equations and the variational principle from which they are 
derived. The novelty of the underlying concepts gives rise to a number of issues and observations. 

PACS numbers: PACS numbers: 98.80.Jk, 98.80.Cq, 95.35.+d, 95.36.+X 



I. EINSTEIN'S UNJUSTIFIED ASSUMPTION 

Albert Einstein, in his 1916 paper Die Grundlage der 
allgemeinen Relativitatstheorie [1] that gave a thorough 
presentation of the theory of gravity he had worked out 
over the preceding decade, made an assumption that does 
not hold up well under scrutiny. Stripped down to its 
barest form the assumption is that inertial mass, and by 
extension energy via E = mc^ , is a source of gravity and 
must therefore be coupled to the gravitational potential 
in the field equations of the general theory of relativity. 
The train of thought that brought him to this conclusion 
is described in §16, where he sought to extend his field 
equations for the vacuum, Rap — ^Rgap = as currently 
formulated, to include the contribution of a continuous 
distribution of gravitating matter of density n, in anal- 
ogy to the extension of the Laplace equation = 
for the newtonian gravitational potential (f) to the Pois- 
son equation V^c/) = Attk^, where k is Newton's gravita- 
tional constant. Einstein referred to fi as the "density of 
matter" , without specifying what was meant by 'matter' 
or its 'density'. Invoking the special theory's identifica- 
tion of "inert mass" with "energy, which finds its com- 
plete mathematical expression in . . . the energy-tensor" , 
he concluded that "we must introduce a corresponding 
energy-tensor of matter " . Further describing this 
energy-tensor as "corresponding to the density /z in Pois- 
son's equation" , he arrived at the extended field equa- 
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tions Rai3 — ^Rgaji = ^§r^Tai3, in which, for a "friction- 
less adiabatic fiuid" of "density" ^, pressure p (a form of 
kinetic energy), and proper four- velocity distribution m"^, 
he took T"^ to be + {p/c'^){u°'u^ - g"^). 

Clearly, Einstein's procedure fails to distinguish be- 
tween the 'active gravitational mass' of matter, which 
measures how much gravity it produces and is the sole 
contributor to the "density of matter" in Poisson's equa- 
tion, and the "inert mass" of matter, which measures how 
much it accelerates in response to forces applied to it, an 
effect entirely different from the production of gravity. 
These two conceptually different masses, along with yet 
a third, all occur in Newton's gravitational equation 

mittB = Fab = -n—^ — , (1) 

in which A/a is the active gravitational mass of a gravi- 
tating body A, mi is the inertial ("inert") mass of a body 
B being acted upon by the gravity of A, and rrip is the 
passive gravitational mass of B, a measure of the strength 
of B's 'sensing' of the gravitational field around A. That 
in suitable units m-i = rup for all bodies is another way 
of saying that all bodies respond with the same accel- 
erations to the same gravitational fields, that, in conse- 
quence, the notion of a 'gravitational force' is irrelevant, 
but the notion of a 'gravitational field' is not. Simple 
thought experiments of Galileo (large stone and smaller 
stone tied together) [2] and Einstein (body suspended 
by rope in elevator) [3] make it clear that bodies do all 
respond alike — an observation now treated as a prin- 
ciple, the (weak) 'principle of equivalence', experimen- 
tally, if somewhat redundantly, well confirmed. That this 
passive-inertial mass has any relation to active gravita- 
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tional mass is not apparent in Eq. (1), where, unlike mj 
and nip, represents a property of A, not of B. But 
Newton's equation for the gravitational action of B on A 
reads 

MiaA = FBA = -«^^. (2) 

Application of Newton's law of action and reaction allows 
the inference that Fab and Fba have the same magni- 
tude, from which follows that ma/mp = Mn/Mp, hence 
that the ratio of active gravitational mass to passive grav- 
itational mass, thus to inertial mass, is the same for all 
bodies. It would seem likely that Einstein relied, cither 
consciously or unconsciously, on this consequence of New- 
ton's laws when he assumed that "inert mass" should 
contribute to the "density of matter" as a source of grav- 
ity in the field equations. 

Newton's law of action and reaction is applicable to 
the bodies A and B only under the condition that grav- 
ity acts at a distance instantaneously, that is, at infinite 
propagation speed. But the general theory of relativity 
that Einstein was expounding is a field theory in which 
gravitational effects propagate at finite speed. Within 
his own theory of gravity there is, therefore, no obvi- 
ous justification for Einstein's assumption that inertial 
mass (and therefore energy) is equivalent to active grav- 
itational mass. This, however, is not to say that there is 
no relation at all between the two kinds of mass. There 
is, for example, the seemingly universal coincidence that 
wherever there is matter made of atoms there are to be 
found both inertial mass and active gravitational mass. 
Indeed, the fact that Newton's theory is a close first ap- 
proximation to Einstein's would argue for some propor- 
tionality between ma and rup for such matter in bulk 
— not, however, for each individual constituent of such 
matter. A 1986 analysis of lunar data concluded that the 
ratio of Wa to mp for aluminum differs from that for iron 
by less than 4 x 10~^^ [4]. An earlier, Cavendish balance 
experiment had put a limit of 5 x 10~^ on the difference 
of these ratios for bromine and fluorine [5]. But those 
results are only for matter in bulk, that is, matter made 
of atoms and molecules. It is entirely possible that elec- 
trons, for example, do not gravitate at all, for no one 
has ever established by direct observation that they do, 
nor is it likely that anyone will. There is in the litera- 
ture an argument that purports to show that if the ratio 
m^/mp is the same for two species of bulk matter, then 
electrons must be generators of gravity [6], but that ar- 
gument can be seen on careful examination to rest on an 
unrecognized, hidden assumption, namely that, in sim- 
plest form, the; gravitational field of a hydrogen atom at 
a distance could be distinguished from that of a neutron 
at the same distance — another assumption no one has 
tested or is likely to test by direct observation. 

Einstein's assumption that energy and inertial mass 
are sources of gravity has survived to the present virtu- 
ally unchallenged. It has generated a number of conse- 
quences that have directed much of the subsequent re- 



search in gravitation theory — indeed, misdirected it if 
his assumption is wrong. Among them are the following: 

• The impossibility, according to Penrose-Hawking 
singularity theorems, of avoiding singularities in 
the geometry of space-time without invoking 'neg- 
ative energy', which is really just energy coupled to 
gravity with polarity opposite to that of the cou- 
pling of matter to gravity. 

• The presumption that the extra, fifth dimension in 
Kaluza-Klein theory must be a spatial dimension 
rather than a dimension of another type. 

• The belief that all the extra dimensions in higher- 
dimensional theories must be spatial, causing the 
expenditure of much effort in mental gymnastics to 
explain why they are not apparent to our senses in 
the way that the familiar three spatial dimensions 
are. 

Denying Einstein's assumption relieves one of the burden 
of these troublesome conclusions and opens the door to 
other, more realistic ones. 



II. NEW, IMPROVED FIELD EQUATIONS 

If Einstein's assumption is to be disallowed, then his 
source tensor for a continuous distribution of gravitating 
matter, T"^ = jiu^u^ + {p/c^){u"'u'^ - g""^), must be 
modified or replaced. One might think to simply drop 
the second term and take T"^ = fiy/^n/^, the cinergy- 
momentum tensor of the matter. This would be incon- 
sistent, for the n in that tensor is the density of inertial- 
passivc mass, which we are now not assuming to be the 
same as active gravitational mass. What to do instead? 

At the same time that Einstein was creating his field 
equations, Hilbert was deriving the field equations for 
(in particular) empty space from the variational principle 
5jR l^l 2 = [7]. This is the most straightforward ex- 
tension to the general relativity setting of the variational 
principle S J |V0p dPx — 0, whose Euler-Lagrange equa- 
tion is equivalent to the Laplace equation V'^cj) = for 
the newtonian potential </>. Modifying that principle to 
S J (|V(/)p -I- 8n Kficj)) (Px = generates the Poisson equa- 
tion V^(/) = AiTKfi. The most straightforward extension 
of this modified principle to general relativity is 

5j{R-^,,)\g\id^ = 0, (3) 

for which the Euler-Lagrange equations are equivalent to 

Ra0 - \Rg<xp = -^i^gaff, (4) 

which makes = —^fiQap, with n now the active 
gravitational mass density, as it should be. Equivalent 
to this equation is Ra/s = ^^IJ-Qap^ the 00 component 



3 



of which reduces in the slowly varying, weak field ap- 
proximation precisely to the Poisson equation. The van- 
ishing of the divergence of the Einstein tensor field on 
the lefthand side of Eq. (4) entails that = Tq,'^.^ = 

— i^in.pga^ + 1^-9 -.p) = —\lJ'.a, hcncc that /X is a con- 
stant. This would seem to be a comedown from the 
equations of motion of the matter distribution implied 
by the vanishing of the divergence of Einstein's T"^ , but 
those equations are unrealistic in that they have the den- 
sity of active mass playing the role that properly belongs 
to the density of inertial mass as the coefficient of the 
four-acceleration of the matter. As will be seen, the new, 
improved field equations to be arrived at will allow ji to 
vary. 

To widen the range of space-time geometries admitted 
by the field equations one can in the usual way add to the 
action integrand of Eq. (3) terms related to such things 
as scalar fields and electromagnetic fields. In particular, 
one can add a cosmological constant term, changing the 
integrand to i? — ^^/x + 2 A and the field equations to 

RafS - ^Rgal^ = -^{fJ' + fi)9oc0, (5) 



where ^^fi = —A. A positive cosmological constant A 
thus appears in this context to be a (mis)representation 
of a negative active mass density /U of a continuous distri- 
bution of gravitationally repulsive matter. An excess of 
this negative density over the positive active mass den- 
sity fi of attractive matter could drive an accelerating 
cosmic expansion, and in doing so would provide a so- 
lution to the vexing 'Cosmological Constant Problem'. 
Leaving aside for the moment the question of how such 
a negative mass density might arise, let us explore the 
consequences of presuming one exists, by studying cos- 
mological solutions of field equations that incorporate a 
positive mass density /i, a negative mass density ft, and 
minimally coupled scalar fields (f> (not the newtonian (p) 
and tl). The variational principle 



S / [R-^{n + fi) + 2(j)-^<j),^-2ilj-'yil),^]\g\U^ = (6) 



combines these elements and generates the field equations 



Ra^ - \Rga0 = ^T^P ■■= -^il^ + ii) gap - 2 {ct>.a<f>.0 - l<l>-^(f>.-y 9^/}) + 2 (V'.a^./3 - ^t^'^.^ 9af3), (7) 

and 

2 {n<l>)cp.„ - 2 (□V)V.a := 2 cp%cp.^ - 2 V^'T^V.a = -^{n + ft). a, (8) 



the latter of which is a derivation from the former, cor- 
responding to the identity T^^ .p =^ 0. In this variational 
principle only the metric components are varied. A 
superficial reason for disallowing variation of (p and tp 
is that varying them would produce the conclusion that 
□0 = and = 0, thus, in light of Eq. (8), the un- 
wanted restriction that /i -h /i is constant. A more fun- 
damental reason for not varying them will be discussed 
later. Notice that the presence of two scalar fields al- 
lows one to disregard by circumvention the question of 
whether the coupling of a scalar field to the space-time 
geometry should have positive polarity or negative polar- 
ity. Indeed, both polarities are consistent with Einstein's 
assumption as well as with its denial, inasmuch as the 'en- 
ergy' of a scalar field x included in x.aX./? - ^X'^X.t 9a/3 
is of a nature entirely different from that of the kinetic 
energy in Einstein's "energy-tensor" , thus should have no 
connection to that tensor. 



and R in centimeters) and dimensionless scalar fields i 
a{t) and ip = P{t) the field equations (7) reduce to 



3^^ = -^(M + M) (9) 



and 



' =-- + + — o^, (10) 



where k = —1, 0, or 1, the uniform curvature of the 
spatial metric ds^. These equations, which are replace- 
ments for the well-studied Friedmann-Lemaitre cosmo- 
logical equations, are equivalent together to 



Attk , , k 6? - , , 



III. COSMIC EVOLUTION EQUATIONS 



and 



For a Robertson- Walker metric (?d1? — l(^{t)ds^ (with 
t in seconds, s dimensionless, c in centimeters/second. 



'^R 



-^(M + M) + ^^. (12) 
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Equations (8) reduce to 



- [a + 3-a]a-^{ (3 + 3-/3] 13 = -—dt{n + fi) 



(13) 

for the time component and da{l-t + fi) = for the space 
components. Equation (13) is equivalent to 

[R%d' - = -^R%{l^ + fi). (14) 

To give substance to the densities /z and p,, let us invoke 
two kinds of gravitating matter: a) primordial matter, 
existent at all times, never changing in amount, with its 
net density ^+p, inversely proportional to the cube of the 
scale factor R; and b) continuously created (or destroyed) 
matter, coming steadily into existence (or passing out of 
existence) at a rate just sufficient to keep its net density 
constant. For primordial matter we have 

R^Uq) 

Mp(i) + Mp(i) = (mp,o +MP,o)-^5^, (15) 

where io denotes the value of t at the present epoch, and 
jjpfi and jlpfi are the present values of the primordial 
densities fxp and p,p. For continuously created matter we 
have 



/«c(i) + /"c(*) = Mc,o + Mc,o, 



(16) 



where /ic,o and p,c,o are the present values of the continu- 
ously created densities /ic and /xc- Under the assumption 
that the net densities fip + ftp and /xc + i^C are additive, 
Eq. (14) turns into 

^dt(R'{d^ - P^)) = ^(/^p,o + /ip,o) i?'(to) at(i?3), 

(17) 

which integrates to 

-[c? - /32) = ^(mp,o + MP.o)-^ + (18) 

where B is the constant of integration, with units cm^. 
At this point it is convenient to set 



:= - 



47r«; 



(AtP,o + A*p,o)-R^(io) 



(units cm) and 



47rK _ 
Ac ■■= 5-(/^c,o + lic,o) 



(units cm ^). Equation (18) becomes 



B - ApR^ 
RP 



(19) 



(20) 



(21) 



and substitution from Eqs. (19), (20), and (21) into 
Eqs. (11) and (12) produces 

^ = ^_A,^__^ (22) 

C2 i?2 3^6 3 ^2 + 3^3 3^6 ) 



and 



1 R _ P2{R) _ Ac Ap 2B 
'^R ~ 3E6 ~ ~3~ ~ 3E3 3E6' 



(23) 



where 



and 



Pi (R) := AcR^ - 3kR'^ + 2ApR^ - B (24) 



P2{R) := AcR^ - ApR^ + 2B. 



(25) 



According to Eqs. (19) and (20), positive values for 
Ap and Ac equate to negative net densities for primor- 
dial matter and continuously created matter, signifying 
that on balance each kind gravitationally repels all other 
matter more strongly than it attracts other matter. This 
excess of repulsion over attraction will be presumed to 
exist for both types. With thus Ap > and Ac > 
assumed, along with an added assumption that B > 0, 
several properties of the scale factor i? as a solution of 
these equations can be inferred rather easily, to wit: 

• For each of fc = —1, 0, 1 there is a positive num- 
ber i?niin (the least positive root of the polynomial 
Pi(i?), where R = 0) that Eq. (22) does not allow 
R to go below. (See Fig. 1.) This sets a posi- 
tive lower bound on the compression of space and 
thereby rules out the development of a 'big bang' 
singularity. 

• For k = 1 and some choices of ^c, Ap, and B 
there are in addition to i?min other values i?max 
and i?min* of R at which Pi, and consequently R, 
vanish. 

• Wherever Pi(i?) < 0, P2(i?) > R^(AcR^ - 2kR + 
Ap) = ipP{(i?). Because, with one exception 
(the graph labeled b in Fig. 1), Pi'(Pmin) > 0, 
also P2(i?min) > 0, and therefore R is positive 
when R = Rmin- This tells that in place of a 
'big bang' there is a 'bounce' off a state of max- 
imum compression at the time when R = Pmin- 
The lone exception has k — 1, P{(i?min) = 0, and 
R{t) = Pmin = Pmax for all time, modeling a static, 
spherical universe of radius Pmin- 

• With the time t — chosen so that i?(0) — Rmin, 
R{t) is symmetric about t = 0, the evolution of 
the universe after time zero thus being mirrored in 
reverse as a devolution of the universe before time 
zero. 
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Pi{R) 
= AqR" - 3kR* + 
{Ac,Ap,B 




fPi(i?) (fc = l) 




FIG. 1: Graphs of Pi{R) for fc = 0, and 1, and generic 
positive values of the parameters Ac,Ap, and B. Values of 
R for which Pi{R) < are excluded from the range of R{t) 
by Eq. (22). 



when k = 1, most notably: (Fig. 1, graph b) static 
spherical universe as noted above; (Fig. 1, graph c) 
spherical universe oscillating for all time between 
minimum radius i^min and maximum radius i?max- 

The polynomial P2{R) is quadratic in R^, with dis- 
criminant 8 Ac S. If B > Al/%Ac, thenP2(i?) 
has no positive root, so i? > at all times, thus the 
universal expansion described by R{t) for times af- 
ter i = is always accelerating. If i? = Ap/SAc, 
then P2{R) has a single positive root, so the ex- 
pansion is, except for a momentary pause, accel- 
erating at all positive times, unless, as in Fig. 1, 
graph b, R{t) = i?min = i?max or, as in Fig. 1, 

graph d, R{t) = i?min* ^max 

then P2(R) has two positive roots. 



liB < A^/^Ac, 



Rc 



which marks a transition from the initial accelerat- 
ing expansion associated with the bounce at t = 
to an interval of decelerating expansion, and 



Ra. '■— 



Ap ^A^-8AgB \ 
2Ac 2 Ac j 



(27) 



which marks a return to accelerating expansion (if 
the expansion is destined to continue forever, as 
opposed to some of the other possible behaviors 
when k = 1). 



Whereas for fc = — 1 and k = (Fig. 1), and k = 1 
(Fig. 1, graph a), the universe expands from i?min 
to infinity as t — )• ±oo, other behaviors are possible 



For the infinitely expanding models the 'Hubble pa- 
rameter' H and the 'acceleration parameter' Q be- 
have asymptotically as follows: 



• 1 2 Ac 3kR^-2At>ll + B Ac /from below if fc > ol 
^ ■= ?^ = ^ 3^^ ' -Y from above if fc < o h' ^ ^''^ 



and 



R/R ^_ 2AcR^-ApR^ + 2B 



2 fcil* - ApR^ + B _ jfrom above if fc > O] 



' + ^ ™ ^ ' from below if fc < ^ ^ ^''^ 



Either of these entails that R{t) Cg^V^c/^''*, • For fc = -1 or 0, and for some cases of fc = 1, 

for some constant C,ast^ ±oo. has a maximum value ilmax a.t R = Rh^^^, where 
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FIG. 2: Graphs of H{R) and Q{R) for generic positive val- 
ues of the parameters Ac, Ap, and B, showing early stage 
inflation foll owed b y a decline in H to its asymptotic limit 
H{oo) = ca/Ac/3, and (for B < Al/8Ac) transitions of Q 
from inflationary acceleration to deceleration at Rd and back 
to acceleration at Ra. to the asymptotic limit Q(oo) = 1. 

dH/dR = c \kR^- ApB?+B)/H R'^ = 0. Iffc = , 
Rh^.. = {^Ja; and ffniax = c v'^c/3 + A^/3B. 
As seen in Fig. 2, H rises sharply from at i?min to 
if max at Rh^^^, then reverses and tails oflF asymp- 
totically to c ^/ Ac/3. One can show that i?inin ~ 
^/B/2A;, Rh^^^ - {^BjA^, Rd - V^BjA^, and 
i^max ~ cAp/^/SB as — >■ with Ac and Ap 
fixed. Thus, as B — ^ with Ac and Ap fixed, 
^min, ^/f,u„x' and i?d are squeezed together closer 
and clos(^r to 0, and iJmax grows without bound. 
This clearly is a recipe for an explosive; postbouncc; 
inflation ending with a 'graceful' exit initiated by 
the onset of deceleration when R reaches R^- 

• As Fig. 2 exhibits, for k = —1, fc = 0, and k = 1, 
graph a, the acceleration parameter Q, desccniding 
from 00 at i?min, passes through 1 at Rh^^^, bot- 
toms out with a minimum value Qmin at J?Q„i„, 
where dQ/dR = 0, then creeps slowly back to 1 as 
i? — > oo (with a late pass through 1 if fc = 1). When 

k = 0, Q^in = -i + ^^JAcB/{Ai,+AcB) at 



FIG. 3: Graphs of H{R) and R{t) for k = 1 and generic 
positive values of the parameters Ac, Ap, and B associated 
with graph c of Fig. 1, showing repetitive, identical periods 
of expansion and contraction, each beginning with a stage 
of rapid inflation from a bounce at i? = Rmin, which is fol- 
lowed by a less rapid expansion to -R = J?max, then a mirror- 
image contraction to an ending stage of rapid deflation into 
the next bounce at 7? = 7?min. The functions are related by 
R{t)/R{t) =: H{R{t)). 

1 

^Q„i„ = {b/Ap + ^B-^jA^ + B/Ac) ' , which, 

with Ac and Ap fixed, is asymptotic to sjBjAc 
as S 0. 

• If fc = 1 and Ac , Ap , and B are such as to generate 
graph c of Fig. 1, then for the section bounded by 
-Rmin and i?max the behavior of R is cyclical, as 
Fig. 3 portrays graphically. 



IV. COSMOLOGICAL SOLUTIONS 
A. Flat open universe (fc = 0) 

When fc = 0, so that space is perfectly flat, fi^in — 

—ApjAc + ■\/Ap + AcB/Ac and it is straightforward to 
integrate Eqs. (22) and (23) (which is equivalent to a 
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For times t w 0, Eq. (31) yields the approximation 



td6 I 1 
^min ' ^ 



3Ac i?;li„ + 



2 i?Ln 

which produces the further approximation that 



(33) 



R{t) « R{t) := iJxnin + 



c 



Rn 



+ 



Ar 



1 



AcR, 



2 

min 



(34) 

thus that R{t) is approximately parabolic, opening up- 
ward from a vertex at .Rminj with an inflationary 'steep- 
ness' factor that grows as as -Rmin 0. 
Equivalent to Eq. (33) is 



R^{t) - R 



3 

min 



3Ar 



R 



3 

min 



(35) 




FIG. 4: Graph of R{t) for fc = and generic positive values 
of Ac, Ap, and B with B < Ap/SAc, showing early accel- 
eration from t = Q to t = td, deceleration from to tg,, 
and ultimately exponential acceleration after ta.; truncated 
graph of R{t) with approximators R{t) and^(t) as identified 
in Eqs. (34) and (36). 



readily integrable equation linear in i?^), with the result 
that 



R'it) 



Ap 

Rmin + -J^ 1 cosh 



At 



SAcct]-^. (31) 



One can show that if B < Ap/8Ac, so that R^ and f?a 
exist, then i?min < Rd ^ Ra- The postbounce times 
and ta. at which R = R^ and i? = i?a are the positive 
solutions of 



cosh I W 3 Ac c t 



^^in + ^PMc ■ 



(32) 



Figure 4 displays the evolution of R{t) for generic values 
of Ac, Ap, and B with B < Ap/8Ac, showing the ac- 
celeration during the interval from t = to t = td, the 
deceleration during the interval from to t^^, and the 
ultimately exponential acceleration after t^,. 



After R{t) has inflated to the point that Raiin/R{t) ^ 1, 
{R^{t) - i?^in)^ wfll difl^er but little from R{t), and then, 
for so long as t remains small enough for the approximate 
Eq. (33) to hold, the approximation 



R{t) f^R{t) := 



3Ac 



^P 
DO _|_ 

^ Ac 



[ctf/^ (36) 



will be valid. Figure 4 shows the approximations to R{t). 



B. Nonflat open {k = —1) and closed (fc = 1) 
universes 

When A: = —1 or 1, elementary symbolic solutions of 
Eqs. (22) and (23) are unavailable, so one must resort 
to numerical solution of Eq. (23). The integration pro- 
ceeds most smoothly when carried out for the equivalent 
equation 



rV o 2 A rr ^C^fc i Sc^Ap 



d2 

-'^min 



min 



(37) 



where U := and S := R/Rmin, with initial condi- 
tions U{0) = 1 and U(0) = 0, obtained from S{0) = 

i?(0)/i?min = Rmin/Rmin = 1 and S{0) = R{0) / R^in = 

0. (Note that this equation and the initial conditions 
apply equally well when fc = 0.) 

To enable the integration, numerical values must be 
chosen for the parameters other than k that appear in 
Eq. (37), namely, Ac, Ap, and B (implicitly through 
-Rmin)- For this determination let us impose the require- 
ment that the solution(s) provide good flts to a Hubble 
diagram that incorporates the data from recent super- 
novae observations. Toward that end a number of defini- 
tions are in order: 

• to, the time of the present epoch, i. e., the time 
elapsed since the bounce at f = 0, in seconds. 



Ho 

rameter 



H{to), the present value of the Hubble pa- 
3r. 
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• Rq := c/Hq, usually referred to as the 'Hubble ra- 
dius' and as the 'radius of the visible universe'. 

• z := R{to)/R — 1, the rcdshift function. 

• Zd := R{to)/Rd — l and z^, := R{to)/R^ — l, the red- 
shifts associated with the transitions from acceler- 
ation to deceleration and from deceleration back to 

acceleration. 

• R := R/ R{to) — 1/(1 + z), a conventional notation. 

• -Rmin '■= -Rmin/-R(io)i a^d -R{d,a} •= -K{d,a}/-R(^o)- 

• A := c/HoR{to) = Ro/R{to), a useful parameter. 



which yields, as roots of ii = 0, 



and 



Rd ■= 



Rn 



2^0 



2nc 



2 fir 



2^0 



(46) 



(47) 



and 



Ok 

Op 

Ob 



H^R^ito) 



SH^R^itoY 



2HoX^ 



-B 



With these definitions Eq. (22) is equivalent to 



1 



1 



= Hnnc + fik— + Op^ + Ob- 

which upon evaluation at to yields 

Oc + 0k + Op + Ob = 1. 
On the other hand, Eq. (23) is equivalent to 

1 



C 



— — — 2 — — 

2 "^R^ i?6 



2 QcR^ - flpR^ - 4 Ob 
2¥ ' 



(38) 
(39) 
(40) 

(41) 

(42) 
(43) 

(44) 
(45) 



(obtainable also from Eqs. (26) and (27) by substitution 
from Eqs. (38), (40), and (41)). 

To produce a curve to match against the data points of 
a Hubble diagram one needs the following standard for- 
mula for the luminosity distance Di, of a photon-emitting 
astronomical object at redshift z: 



Dj^{z) = {l + z)R{to)ri,\ 



R{to) Jo H{u) 



du 



(48) 



Here 



^k(p) := 



sinh/9 if k = — 1, 
p if fc = 0, 

sinp if A; = 1, 



(49) 



and 



H{z):=Ho[nc+n^{i + zf+np{i + zf+nB{i + zf\^. 

(50) 

In terms of the parameter A the luminosity distance for- 
mula reads 



D^{z) = {l + z)—r^[X 



du 



[Oc + Ok(l + uY + Op(l + uf + Ob(1 + w)6] 2 



(51) 



From Dl one constructs the distance modulus /i, the 
difference between the apparent magnitude and the ab- 
solute magnitude of the object in question at redshift z\ 



this reduces in the usual way to 
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Prom R = l/(l + z) follows -z/{l + z) = R/R = H{z), 
which upon integration yields the following formula for 



the time T(z) elapsed between an event occurring at red- 
shift z and the present time to (when z = 0): 



T{z) 



{l + u)H{u) 



■ du — — — 



du 



Ho Jo (l + u)[f2c+f^k(l+w)2 + Op(l + u)3 + OB(l + w)<']5 
I 



(53) 



The distance modulus fi of Eq. (52) is the function to 
be fit to the relative magnitudes data obtained from ob- 
servations of supernovae. The data to be used are those 
from the 182-member sample described in Sec. 3 of [8], 
provided at http://braeburn.pha. jh.u.edu/~ariess/ 
R06/sn_seiiiiple. The requisite calculations are carried 
out (using the indicated Mathematica functions) in two 
stages, the first of which comprises the following steps: 

1. On first run, set Hq = 74.2 km s~"'^/Mpc, the recent 
estimate of Hq obtained by Riess, et al. [9]. On 
subsequent runs set Ho = 70.6 km s~^/Mpc and 
Ho = 77.8 km s'VMpc, the extremes of the 4.8% 
confidence interval predicted in [9]. 

2. Set B = 0, thus Ob = 0. This makes i?d = i?min = 
0, so will produce a 'big bang' model with = oo 
and with no inflation (cf. Eq. (23), in which the 
last term, now absent, is the one that makes R 
large when i? « 0). 

3. For fc = 0, using NIntegrate as the integrator for 
Di,{z) in Eq. (51), apply NonlinearModelFit to 



J 



IJ.{z) to obtain, subject to the constraint of Eq. (43), 
Chi-square best-fit values for ilp, fie, and A. Do 
the same for fc = — 1 and fc = 1. (Note that when 
fc = 0, A disappears from Eq. (51), so is irrelevant 
to the fitting.) 

4. From the best-fit values thus obtained, compute Rg, 
from Eq. (47), then z^, from R^, = 1/(1 -|- Zg,). 

5. From Eq. (53), by applying NDSolve to the initial- 
value problem dT/dz = 1/(1 + z)H{z), T(0) = 0, 
compute T{z) and then the implied time elapsed 
since the big bang as lim T(z). 

6. Compute the Chi-square statistic 



182 . ^2 

^2 ^ / /i(Zj) - fJ.t 



Table I presents the results of these computations. 



(54) 



TABLE I: Best-fit parameters (A not shown) for 'big bang' model: B set to 0, thus fls = 0, so Rd = Rmin = and Zd = oo. 



Inputs 


Outputs 


Ho (km s'VMpc) 


k 




fie 








lim T{z) (Gyr) 

z—^oo 






-1 


0.419 


0.581 


2.986 X 10"^° 


0.712 


0.405 


12.14 


163.075 


70.6 





0.419 


0.581 





0.712 


0.405 


12.14 


163.075 




1 


0.419 


0.581 


-2.379 X 10"^* 


0.712 


0.405 


12.14 


163.075 




-1 


0.295 


0.705 


2.796 X 10"^^ 


0.594 


0.684 


12.77 


160.829 


74.2 





0.295 


0.705 





0.594 


0.684 


12.77 


160.829 




1 


0.295 


0.705 


-2.796 X 10"^^ 


0.594 


0.684 


12.77 


160.829 




-1 


0.198 


0.802 


1.262 X 10"^^ 


0.498 


1.008 


13.57 


183.225 


77.8 





0.198 


0.802 





0.498 


1.008 


13.57 


183.225 




1 


0.198 


0.802 


-1.277 X 10"^^ 


0.498 


1.008 


13.57 


183.225 
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To prepare for the second stage observe that Eqs. (46) 
and (47) imply that 

np = 2{Rl + Rl)nG and I^b = -5 ^^c- (55) 

Equivalent to Eq. (43) is ilc+f^p + f^B = 1 — f^k, substitu- 
tion into which from Eqs. (55), then back into Eqs. (55), 
yields 



Of 



l + 2{Rl + Rl) 



2 -"-d -"-a 



2(i?3 



^D(l-f^k) 



l + 2{Rl + Rl) 



1 p3 p3 ' 

2 ^d^a 



and 



Or = 



-iR^Rlii-n^) 



1 + 2{RI + Rl) 



2 -"-d -"-a 



(56) 



(57) 



(58) 



where, from Eq. (39), Ok = — A^fc. 

Because the supernovae data extend into the past no 
further than z = 1.755, the best-fit parameters Op, Oc, 
and A will change by negligible amounts when Ob changes 
from to a small nonzero number. For this reason a dif- 
ferent strategy is called for in the second stage, as de- 
tailed in the following steps: 

7. Set A and Za to the best-fit values found in the first 
stage (for k = 0, choose A = 1). Assign zd a large 
value. 



J 



8. From R = 1/ (l+z), compute R^ and i?a (reversing, 
in the latter case, the computation of Za. from i?a 



in the first stage). 



9. Using Eqs. (56), (57), (58), and Ok = -A^A;, com- 
pute Oc, Ok, Op, and Ob, and in turn Dl{z) and 

H{z) from Eqs. (51) and (52), applying NDSolve 

* 

to the initial- value problem dl/dz = Ho/H{z), 
/(O) = to compute the integral in Eq. (51). 

10. Compute the Chi-square statistic and compare it 
to the corresponding number from the first stage. 
If it is appreciably larger, make Zd larger, thus Rd 
smaller and Ob closer to 0. 

11. When a good match is found in Step 10, recover 
Ac, Ap, and B from Eqs. (38), (40), and (41), 
then iJniin from Pi(i?min) = (cf. Eq. (24)) us- 
ing NSolve. 

With Ac, Ap, B and i?min in hand, proceed to the 
solution by NDSolvc of Eq. (37) for Uif) and S{f) = 
[U{t)]^, with C/(0) = 1 and C/(0) = 0, to obtain a fuh 
history of the model universe. As a check use FindRoot 
to get from S{to)/S{to) = -ff(to) — -^o and compare it 
to T(zmin) from Eq. (53) (see Step 5), where 1 -|- Zmin = 
R{to)/Rmin = Ro/XRmin- They must agree. 

One can find iioo, the time at which R has doubled 
in size one hundred times from its minimum size at the 
bounce, by solving the equation log2 S{t) = 100 using 
FindRoot. If a different value for iioo is wanted, change 
Zd and repeat Steps 7-11 and the integration of Eq. (37). 

Table II presents the results of these computations. 



TABLE II: Best-fit parameters for 'bounce' model: B > 0, thus SIb < and -Rmin > 0. (The columns for Q,p, Q,c, and i\ are 
identical to those of Table I, and are included by reference.) 



Inputs 


Outputs 


Ho (km s-VMpc) 


k 


A 






Zd 


Ob 




iloo (s) 


to (Gyr) 






-1 


1.728 X 10" 


20 


0.405 


1 X lO^'' 


-1.047 X 10" 


196 


1.02 X 10"^^ 


12.14 


163.075 


70.6 





1 




0.405 


1 X lO^-"' 


-1.047 X 10" 


196 


1.02 X 10"'^'''' 


12.14 


163.075 




1 


1.543 X IQ- 


-9 


0.405 


1 X lO''^ 


-1.047 X 10" 


196 


1.02 X 10"^^ 


12.14 


163.075 




-1 


1.672 X 10" 


11 


0.684 


1 X lO"'^ 


-7.378 X 10" 


197 


1.15 X 10"^^ 


12.77 


160.829 


74.2 





1 




0.684 


1 X lO^-"' 


-7.378 X 10" 


197 


1.15 X 10"^^ 


12.77 


160.829 




1 


1.672 X 10" 


11 


0.684 


1 X lO^'' 


-7.378 X 10" 


197 


1.15 X 10"^^ 


12.77 


160.829 




-1 


3.552 X 10" 


10 


1.008 


1 X lO"'' 


-4.951 X 10" 


197 


1.34 X lO"^'^ 


13.57 


183.225 


77.8 





1 




1.008 


1 X 10"-' 


-4.951 X 10" 


197 


1.34 X 10"^^ 


13.57 


183.225 




1 


3.573 X 10" 


10 


1.008 


1 X lO^'^ 


-4.951 X 10" 


197 


1.34 X 10"^^ 


13.57 


183.225 
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FIG. 5: Best fit models of n(z) for fe = (parameters shown in Tables II and I), plotted against the 182 SNe la data points 
they are fitted to. The corresponding plots for k = —1 and k = 1 are visually indistinguishable from these. 



Figure 5 shows the best-fit models of ^{z) for k — 0, 
plotted against the 182 SNe la data points they are fit- 
ted to. The parameters are those of Table II. The corre- 
sponding plots for k = —1 and k = 1 are visually indis- 
tinguishable from the ones shown. 

Figure 6 exhibits a complete postbounce history of 
the best-fit model for Hq — 74.2 km s^^/Mpc and 
k = 0. The graph is constructed from R^^^ = 

-Ap/Ac + ^/A^+AcB/Ac, 5* := i?/i?mi„, and the 
equation, equivalent to Eq. (31), 

RHt) = ^^^^^ cosh (VM^ca)-^, (59) 

Ac V / Ac 

with Ap = 5.519 x 10^^ cm, Ac = 1.360 x lO"^^ cm'^, 
and B = 5.348 x 10^'*'* cm'', obtained from the best-fit 
parameters fip, fie, and f^B by use of Eqs. (40), (38), 
and (41). The history comprises a brief interval of rapid 
inflation from the bounce, sliding into deceleration and a 
graceful exit into an era of coasting, followed by a return 
to ultimately exponential acceleration at z = Za = 0.684 
(when t — 6.90 Gyr). Throughout the coasting era 
Qit) ~ Qmin = -0.5 + 7.33 X 10-^^, reflective of the 
approximate proportionality R(t) oc {ct)"^^^ implied by 
Eq. (36). At t = t^, Q{t) = 0.557. Numerical integration 
of Eq. (37) produces the same graph for fc = 0, and for 



fc = ±1 yields graphs visually indistinguishable from it. 



V. DARK MATTER AND DARK 'ENERGY' 

A. Drainholes 

For their success in modeling infiation, deceleration, 
and acceleration the cosmological solutions developed in 
Sec. IV depend on, among other modifications, the inclu- 
sion of the negative density /2 alongside the usual positive 
density fi in the variational principle of Eq. (6) and the 
consequent field equations. To make these models cred- 
ible requires finding a source for the implied repulsive 
gravity. To attribute the repulsive gravity to a mysteri- 
ous dark 'energy' represented by a cosmological constant 
A would be to beg the question and would add nothing to 
our understanding of the cosmos. As it happens there is 
ready to hand a well-aged theoretical discovery that has 
negative, repulsive gravity as an integral feature. In 1973 
I described in considerable detail a model of a gravitating 
particle alternative to the Schwarzschild vacuum solution 
of Einstein's field equations. This space-time manifold, 
which I termed a 'drainhole' with 'ether' fiowing through 
it, was discovered independently at about the same time 
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.S 
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FIG. 6: Postbounce graph of logj S{t) versus log2(l + f ) [t := t s^^) for the best-fit solution with Ho = 74.2 km s~^/Mpc 
and fc = 0. The early stage rapid inflation, after producing 177 doublings of the normalized scale factor 5 in the first 
second after the bounce, gives way smoothly to a long period of uphill, dccclerative 'coasting' (where the graph is nearly 
linear). Initially, the acceleration parameter Q := {R/R)/[R/Rf = {S/S)/{S/Sy is positive and huge ((3(0) — oo). In 
the coasting period Q{t) w Qmin = —0.5 + 7.33 x 10^® , which is reflected in the observation that in that interval logj S{t) 
177+((217-177)/(58-0))log2(l+t) = 177+(20/29) log2(l+f ) « 177+(2/3) log2(l+t), so that S{t) « 2^''''{l+ty 



2/3 



177^2/3 



(5(t) oc f^^"^ Q(i) = —1/2). After the coasting era Q(t) becomes positive again at redshift Za = 0.684 (when t = 6.60 

Gyr and log2(l + t) — 57.6), rising to 0.557 at t = to, and then asymptotically to 1 as i — > oo, the growth rate of S{t) thus 
becoming asymptotically exponential. The corresponding graphs for the best-fit solutions with Ho the same and k = ±1 are 
visually indistinguishable from this one. 



by K. A. Bronnikov and formally some years earlier by O. 
Bergmann and R. Leipnik who rejected it for "physical 
reasons" [10-12]. It has come to be recognized as an early 
(apparently the earliest) example of what is now called 
by some a 'traversable wormhole', and has been analyzed 
from various perspectives by others [13-27]. The met- 
ric is a static, spherically symmetric solution of the field 
equations (7) and (8) with jj. = j2 = and ■0 = 0. (NB 
Rafi and R here are the negatives of those in [10].) It 
has the proper-time forms (in units in which c = 1) 

dr' = [1 - f{p)] dT' - [1 - f{p)]-' dp' - r'ip) dn' 

(60) 

= df -[dp- f{p)dt]'^ -r'{p)dn', (61) 
where i = T- j f{p)[l - p{p)]-^ dp, 

/2(;9) = l-e-(2'"/"'"(''\ (62) 



r{p) = V(p - m)2 + a2 e^™/" 



and 



(j) = a{p) = - 



tan^i 

2 \ a 



(63) 



(64) 



with a := — m?, the parameters m and n satisfying 
< m < n. (The coordinate p used here translates to 
p+TO in [10].) The shapes and linear asymptotes of r and 
P are shown in Fig. 7. Not obvious, but verifiable, is that 
P{p) ~ 2m/ p as p — > oo, which, together with r(p) ^ p 
as /9 — > oo, shows m to correspond to the Schwarzschild 
mass parameter. 

The shared metric of the cross sections of constant t is 
dp'^ + r'{p) dfl"^. The choke point of the drainhole throat 
is the 2-sphere at p = 2m, of areal radius r(2m) = 
j^g(m/n)a(2TO)^ which incrcascs monotonically from n to 
ne as m increases from to n. Thus the size of the throat 
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is determined primarily by n, and only secondarily by to. 
Although the scalar field (p has a nonkinetic 'energy' den- 
sity that contributes to the space-time curvature through 
Tap, this energy has little to do with the strength of grav- 
ity (as determined by to), rather is associated with neg- 
ative spatial curvatures found in the open throat, cur- 
vatures whose negativity mandates the unconventional 
choice for the polarity of the coupling of (f) to the space- 
time geometry, as expressed by a minus sign in X'q,^ in 
Eq. (7). Even when to = 0, so that there is no gravity, 
these negative curvatures persist and the throat remains 
open. As a matter of perspective, it is more insightful 
to consider that the scalar field does not cause (i. e., is 
not a source of) these spatial curvatures, but simply tells 
of their existence and describes their configuration. This 
perspective helps disabuse one of the peculiar notion that 
geometrically unexceptionable space-time manifolds such 
as the drainhole are somehow a product of 'exotic' matter 
just because their Ricci tensors disrespect some 'energy 
condition'. Moreover, it is not a great stretch to sur- 
mise that, whereas the parameter to specifies the active 
gravitational mass of the (nonexotic) drainhole particle, 
the size parameter n specifies in some way its inertial 
rest mass. This speculation is supported by two consid- 
erations: first, as shown in [10], the total energy of the 
scalar field (j) lies in the interval from n/2 to n7r/2, thus is 
essentially proportional to n; second, it would seem likely 
that the bigger the hole, the greater the force needed to 
move it. A 'higgsian' way of expressing this idea is to 
say that the drainhole 'acquires' (inertial) mass from the 
scalar field (j). 

Because r{p) > n > and /^(p) < 1, the drainhole 
space-time manifold is geodesically complete and has 
no one-way event horizon, the throat being therefore 
traversable by test particles and light in both directions. 
The manifold is asymptotic as p — > oo to a Schwarzschild 
manifold with (active gravitational) mass parameter to. 
The flowing 'ether' (a figurative term for a cloud of iner- 
tial observers free-falling geodesically from rest at p = oo) 
has radial velocity f{p) (taken as the negative square root 
of P{p)) and radial acceleration (/^/2)'(p), which com- 
putes to ~m/r'^{p) and therefore is strongest at p = 2 m, 
where r assumes its minimum value. Because the ra- 
dial acceleration is everywhere aimed in the direction 
of decreasing p, the drainhole attracts test particles on 
the high, front side, where p > 2 m, and repels them 
on the low, back side, where p < 2 to. Moreover, the 
manifold is asymptotic as p ^ — oo to a Schwarzschild 
manifold with mass parameter to = — TOe™''/", so the 
drainhole repels test particles more strongly on the low 
side than it attracts them on the high side, in the ra- 
tio —to/to = e^'^Nn-^-m.'\ r^Yie drainhole is a kind of 
natural accelerator of the 'gravitational ether', drawing 
it in on the high side and expelling it more forcefully on 
the low side. To replace the somewhat disreputable term 
'ether' with something more acceptable in polite society 
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FIG. 7: Graphs of r(p) and /^(p) for gene ric value s of the 
parameters m and n {0 < m < n and a := — m^). 

one can imagine that it is space itself that is flowing into 
and through the drainhole. This substitution, which ac- 
cords with Einstein's insight that the concepts of space 
and of a gravitational ether are essentially interchange- 
able [28], should cause no alarm, for the very notion of 
an expanding universe already ascribes to space the req- 
uisite plasticity. 

The discovery of the drainhole manifolds arose, in my 
case, from a search for a model for gravitating particles 
that, unlike a Schwarzschild space-time manifold, would 
have no singularity. Geodesic completeness and absence 
of event horizons resulted naturally from that require- 
ment and the invoking of a minimally coupled scalar 
field to weaken the field equations. As shown in [10], 
a drainhole possesses all the geodesic properties that a 
Schwarzschild blackhole possesses other than those that 
depend on the existence of its horizon and its singular- 
ity, having evicted the horizon and replaced the singu- 
larity with a topological passageway to another region of 
space. Drainholes are able, therefore, to reproduce all 
the externally discernible aspects of physical blackholes 
(if such things exist) that Schwarzschild blackholes repro- 
duce. That their back sides have never been knowingly 
observed (but in principle could be) is no more troubling 
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than the impossiblity of directly observing the back sides 
of Schwarzschild horizons. For these reasons drainholes 
are more satisfactory than Schwarzschild blackholes as 
mathematical models of centers of gravitational attrac- 
tion. Moreover, there is little reason to doubt that ro- 
tating drainhole manifolds analogous to the Kcirr rotat- 
ing blackhole manifolds exist and will prove to be better 
models than the Kerr manifolds. (A relatively recent so- 
lution of the field equations (7) and (8) perhaps describes 
such a manifold [26].) 

B. Dark matter and dark 'energy' from drainholes 

A physical center of attractive gravity modeled by a 
drainhole would qualify to be called a 'darkhole', inas- 
much as (as shown in [10]) it would capture photons that 
venture too close, but, unlike a blackhole, must even- 
tually release them, either back to the attracting high 
side whence they came or down through the drainhole 
and out into the repelling low side. Thus one can imag- 
ine that at galactic centers will be found not supermas- 
sivc; blackholes, but supermassivc; darkholes. At the cen- 
ter of our galaxy, for example, instead of a blackhole 
of Schwarzschild mass m « 4 x IO^Mq « 1.18 x lO'^ 
km (in geometric units), and of horizon area 47r(2 m)^ » 
4.39 X 10^^ km^, there might reside a drainhole the area 
47rr^ (2 m) of whose throat would lie between 47rn^ and 
47r(ne)^, with n constrained only by having to exceed 
m. As measured by the bounds n and ne on the areal 
radius r{2m), and the corresponding bounds —m/ATTV? 
and —m/An{nef' on the maximum radial acceleration 
(/^/2)'(2 to) — — m/r^(2TO), such an object might be 
large but relatively benign (n^O,TO/n<Cl), large and 
not so benign (n » 0, m/n ~ 1), small and malignant 
(n ~ 0, m/n « 1), small and not so malignant (n ~ 0, 
m/n <C 1), or something in between. 

There is, however, more to the story. A central tenet 
of the general theory of relativity is that every elemen- 
tary object that 'has gravity' is a manifestation of a local 
departure of the geometry of space-time from flatness. 
If such an object has other properties ascribed to it by 
quantum theory, these must be additional to the underly- 
ing geometric structure. I therefore propose the hypoth- 
esis that every such elementary gravitating object is at 
its core an actual physical drainhole — these objects to 
include not only elementary constituents of visible mat- 
ter such as protons and neutrons, or, more elementarily, 
quarks, but also the unseen particles of 'dark matter' 
whose existence is at present only inferential. Moreover, 
I will consider that visible matter and primordial matter 
are, gravitationally, one and the same, and that the dark 
matter is the continuosly created matter. 

The pure, isolated drainhole described by Eqs. (60-64) 
is an 'Einstein-Rosen bridge' between two otherwise dis- 
joint 'subuniverses', each of which, if evolving, would by 
itself consume for its description all the resources of a 
Robertson-Walker metric. Nonisolated drainholes could 



presumably exist not only as 'bridges' between our sub- 
universe and another (or multiple others), but also as 
'tunnels' from one place in our subuniverse to another, 
possibly quite distant from the first by every nontunnel 
route. Both types could contribute to the negative mass 
density jl as well as to /z, the bridge drainholes contribut- 
ing to jl if their gravitationally repulsive back sides re- 
side in our subuniverse, the tunnel drainholes contribut- 
ing to /i by way of their gravitationally repulsive exit 
portals. Tunnel drainholes are easy enough to visualize 
in abundance as topological holes into which space dis- 
appears, only to reappear elsewhere in our subuniverse, 
in analogy with rivers that go underground and surface 
somewhere downstream. Having both portals located in 
our subuniverse, tunnel drainholes would have proper- 
ties we could fully take account of. Bridge drainholes, 
on the other hand, with only one; side; in our subuni- 
verse would have properties dependent in part on cir- 
cumstances in other subuniverses, circumstances beyond 
our ken. I shall, therefore, assume that no bridge drain- 
holes contribute to jl (or to /z) , that the only contributors 
are tunnel drainholes, consequently that our subuniverse, 
even if part of a 'multiverse', is self-contained. 

Lacking for the present a full mathematical description 
of these tunnel drainholes, let us nevertheless proceed as 
if they exist and are characterized by parameters to and n 
related as in an isolated bridge drainhole. Every particle 
of gravitating matter, whether in the 'primordial' (P) 
category or the 'continuously created' (C) category, is 
then at its core one of these tunnels, and it becomes a 
question of relating m and n to the densities /Up,o, Mp,o, 
Hcfi, and jlcfi- 

The fitting of the cosmological solutions of Sec. IV to 
the SNe la data determines, by way of Eqs. (19), (20), 
(38), and (40), and Table II, the net densities /Up,o + A*p,o 
and A*c,o + Mc,0; but tells nothing about the mix between 
the positive densities and the negative. To 'fix the mix', 
even if somewhat arbitrarily and only for P-densities, a 
procedure comes to mind that at first glance; appears to 
replace a negative net density with a positive density, but 
on closer inspection is seen not do so. Equation (42) can 
be recast in the form 

tt2 Sttk, _ , 1 

H = y-4c - — (mp,o + MP,o)-^ 

+ ifo^(ai^+OBi^). (65) 

The analogous equation for the 'concordance' (ACDM) 
model is 



+ J?o2(a-i^ + f}R-ij), (66) 

where the subscripts M and R refer to matter and radia- 
tion, respectively. As the only kind of gravitating matter 
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contemplated by the concordance model is effectively pri- 
mordial matter (existent from the epoch of baryon gen- 
esis or earlier, treated in Eq. (66) as if existent from the 
big bang onward), the density /xm can be identified with 
the density fJ.p,o, and then the corresponding primordial 
matter terms in the two equations will match exactly if 
-(/ip,o + Ap,o) = Mm = Mp,o, thus if /<p,o/Aip,o = -2. 

To produce P-densities in the ratio flp / with tunnel 
particles of uniform masses mp and sizes np requires that 
flp /hp = mp/mp = -e'"i"'/V»l-'"i>, which entails that 
p,p/Hp is constant and 



mp 
np 



[ln(-^p//ip)] 
[ln(-Aip/Mp)]' ^ 



r2 • 



(67) 



If indeed /ip,o/MP,o = —2, then ^p/fip = —2, and 
mp/np = In2/^(ln2)2 -K7r2 ^ o.22. Under the assump- 
tion that np, or better, the throat radius r(2mp), is a 
proportional measure of the inertial mass of the primor- 
dial particles, this fixing of mp/np would provide a basis 
outside the context of newtonian physics for suspecting 
the existence of a uniform ratio of active gravitational 
mass to passivc-incrtial mass of those particles. If, as 
I propose, the primordial particles are taken to be the 
gravitational cores of the protons and neutrons of the 
present epoch, then perhaps r{2mp) = mp_n, the inertial 
mass of a proton or a neutron. Converted to geometric 
units mp.n = 1.67 x lO'^^ g = 1.24 x 10"^^ so that 
np = r(2TOp)e-('"p/"p)"(2™p) = 1.24 X 10-24 g = 9.21 x 
10-^3 cm and mp = 2.67 x lO'^^ g = 1.99 x lO'^^ ^m. 
Prom Eqs. (19) and (40) one calculates that 



Atp.o + MP,o = — 



3Hi 



Sttk 



Op 



(68) 



which for the best-fit value ftp = 0.295 of Table II (cf. 
Table I) reduces to ^p,o + p-p.o = -3.05 x lO^^*' g/cm^ 
= -2.27 X 10-58 cm/cm^. With /Up,o/mp,o = -2 this 
gives fJLpfi/mp = 11.4/m^ as the present proton-neutron 
{qua primordial particle) number density. 

Comparison of the terms c^Ac/S and c^A/S in 
Eqs. (65) and (66) serves only to provide an identity 
for A, so to 'fix the mix' of positive and negative densi- 
ties in /ic,o + /^C,o requires a strategy different from that 
used for HPfi -\- jlpfi. Note that, unlike the concordance 
model, which combines the gravitational effects of bary- 
onic (visible) matter and cold dark matter in the term 
(87rK/3)(/UM/-R^) in Eq. (66) and thereby predicts from 
observational data a cold dark matter (CDM) to baryonic 
matter (BM) mass ratio of approximately 4:1, presum- 
ably constant in time, the model developed here reserves 
the analogous term in Eq. (65) for the effects of baryonic 
matter alone. That leaves the CDM to BM ratio, speci- 
fied here as ^c/fJ-p, undetermined, and morever allows it 
to be time-dependent, inasmuch as He + flc is constant 
and yUp -h/Up oc R~^. To generate estimates for C- matter 
like those for P-matter, let us set 7(2;) := iJ,c{t)/ Hp{t), 



where t is the epoch of events seen at redshift z, and 
use 7o := 7(0) = Hc{to)/ l^p{to) as a parameter. From 
Eqs. (20) and (38) one calculates that 



(69) 



This, together with Eq. (68) and the stipulation that 
A*p,o/a*p,o = ^2, implies that 



Mc.o 
Mc,o 



1 + 



2 0c 
7of^F 



(70) 



If, as in the case of P-particles, we assume that all the 
tunnel C-particles have the same mass parameter mc and 
same size parameter nc, then, just as for P-particles, we 
have that p,c / Mc is constant and 



mc 
nc 



IH-Hc/Hc)? 
[ln(-/xc/Mc)]2 + 



(71) 



Using this equation in conjunction with Eq. (70) and 
the equation nc = r(2 mc)e-'^'"<^/"'^^"(2"*°\ we can con- 
struct a table of output values of parameters of interest 
for selected input values of 70 and r(2mc). The result 
is Table III, which uses Hq = 74.2 km s~^/Mpc and the 
best-fit values Op = 0.295 and Clc = 0.705 from Table II 
(cf. Table I). 

This discussion has assumed that the mass and size pa- 
rameters of the P- and C-particles remain constant over 
time, which need not be the case. Indeed, it is conceiv- 
able, perhaps even likely, that these particles would share 
in the universal contraction and expansion. A particle 
that did so could be represented by a drainhole model 
that over the course of time from —00 to 00 shrinks from 
infinite size to a point and then grows back to infinite 
size. Such a model is presented in [29] as a solution of 
the same field equations that the static drainholes sat- 
isfy, namely, Eqs. (7) and (8) with fj, = p, = and tp = 0. 
This model has m = 0, thus is devoid of gravity, but it 
could in principle be modified to one that has m ^ 0. 
Allowing the P- and C-particlcs to contract and expand 
with the universe would, of course, necessitate a revision 
of the discussion above and of Table III. 



VI. ISSUES AND OBSERVATIONS 

The conceptual basis of the model of the cosmos de- 
veloped in the preceding sections differs in many respects 
from that of the concordance model. These differences 
and their consequences raise a number of issues, some 
of which have been discussed above to greater or lesser 
extent, others of which have been passed over. In this 
section I shall specify and comment on those that I have 
been able to identify along the way, in approximately 
the order in which they appeared. For ease of reference I 
shall call my model the PCB model, and the concordance 
model the ACDM model. 
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TABLE III: Computed parameter values for continuously created C-matter: nc (attractive mass density), mc (attractive mass 
of a particle), /Zc/mc (ratio of repulsive mass density to attractive mass density), and fic/mc (particle number density). Inputs 
are a selection of assumed values of r(2mc) (inertial mass of a particle) and 70 (ratio at present of attractive mass density 
of C-matter to that of primordial P-mattcr, proposed as dark matter to baryonic matter density ratio). The best-fit values 
Op = 0.295 and Oc = 0.705 are taken from Table II (cf. Table I), with Ho = 74.2 km s"VMpc. 



Inputs 


Outputs 


r(2mc) 


Mc,o 

70 := 

A*P,o 


MC (g/cm ) 


mc (g) 


MC 


A*C /, , 3\ 

mc 




10 


3.05 X 10"^® 


1.72 X 10"^^ 


-1.48 


1.77 X 10^ 


rup-n = 1.67 X IQ-^* g 
= 1.24 X 10"^^ cm 


100 
200 


3.05 X 10"^® 
6.10 X 10"^* 


2.43 X 10"^® 
1.24 X 10"^® 


-1.05 
-1.02 


1.26 X 10* 
4.91 X 10* 




■)UU 


i.:)-) X lU 


■ ).U .L X J.U 


I ("1 1 
- i . U J. 


'"> rv* w 1 n' 
■ ).\)') X lU 




10 


3.05 X 10"^^ 


1.72 X 10"^^ 


-1.48 


1.77 


10^ X mp-n 


100 
200 


3.05 X 10"^® 
6.10 X 10"^® 


2.43 X 10"^" 
1.24 X 10"^" 


-1.05 
-1.02 


1.26 X 10^ 
4.91 X 10^ 




500 


1.53 X lO"^'' 


5.04 X 10"^^ 


-1.01 


3.03 X 10^ 




10 


3.05 X 10"^^ 


1.72 X 10"^^ 


-1.48 


1.77 X 10"® 


10^° X mp.„ 


100 

200 


3.05 X 10"^® 
6.10 X 10"^® 


2.43 X 10"^® 
1.24 X 10"^® 


-1.05 

-1.02 


1.26 X 10"^ 
4.91 X 10"^ 




500 


1.53 X 10"^^ 


5.04 X 10"^^ 


-1.01 


3.03 X 10"^ 


mpianck = 2.18 X lO"'^ g 

= 1.62 X 10"^^ cm 
= 1.31 X 10^^ X mp-n 


10 

100 
200 
500 


3.05 X 10"^^ 
3.05 X 10"^* 
6.10 X 10-2® 
1.53 X lO-^'' 


2.24 X 10"^ 
3.16 X 10"'^ 
1.62 X IQ-'^ 
6.56 X 10"® 


-1.48 
-1.05 
-1.02 
-1.01 


1.36 X lO"^'' 
9.65 X 10"^" 
3.77 X 10-1^ 
2.33 X 10"" 



A. Gravity and passive-inertial mciss 



As noted at the beginning, Newton's law of action and 
reaction allows the inference that the ratio rna/mp.i of 
active gravitational mass to passive-inertial mass is the 
same for all 'bodies', that is, for all matter in bulk. Be- 
cause Newton's theory describes gravity as well as it does, 
this inference must be at least approximately correct. 
Even if exact, however, it is only a statement about a nu- 
merical ratio of quantities that conceptually have nothing 
to do with one another: the generation of a gravitational 
field on the one hand, the resistance to being acceler- 
ated by a field of any sort on the other. Tb try to force a 
marriage between these concepts is akin to the proverbial 
task of trying to push a square peg into a round hole. If 
one were to treat this unlikely union as legitimate, then 



one would have to say that a gravitating body flying by 
at speed v would have its active gravitational mass and 
therefore its gravity increased by a factor l/-\/l — v'^/c^ 
because that's what happens to its inertial mass. But 
the 'increase of inertial mass' by that factor has nothing 
to do with an innate property of the body such as the 
strength of its gravity, rather is simply a quantification of 
the fact that the closer any body's speed is to the speed 
of light, the harder it is to accelerate that body with a 
field that propagates at the speed of light. (The analogy 
of a person on foot pushing a dead-battery car to get it 
started comes to mind: if the car can't be started at a 
speed less than the pusher's maximum running speed, it 
won't be started.) And let us note in passing that when 
an electrically charged particle's electromagnetic field is 
observed to depend on the relative velocity of the ob- 
server and the particle, this is attributed to the fact that 
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the field propagates at finite speed, not to a change in 
the particle's electric charge. 

If passive-inertial mass is not a source of gravity, then 
as said in Sec. II Einstein's 'energy-tensor of matter' 
T'^P := ij,u°'uP + {p/c^){u°uP - .g"'^) has no role to play 
in the field equations of gravity, so one loses the spu- 
riously attractive feature that the (Einstein) field equa- 
tions imply that T'q'^.^ = 0, interpreted as the equations 
of motion of the gravitational source matter. This, how- 
ever, does not keep one from asserting that Tq,'^.^ = 0; it 
only requires that one look elsewhere for its justification. 
But in making that assertion one must recognize that, 
for Ta^.^ = to be so interpreted, /U must be the density 
of inertial-passive mass, not of active gravitational mass. 

All too frequently there appear in the popular press ac- 
counts of astronomers 'weighing' the universe or 'weigh- 
ing' a galaxy. One can object, of course, to the pre- 
sumably intentional blurring of the distinction between 
inertial mass and weight, but a larger objection is that 
such accounts leave the reader with the false impression 
that determining an entity's active gravitational mass is 
known to be equivalent to determining its passive-inertial 
mass (and therefore its weight). Only if one can prove 
the existence of a universal proportionality between the 
two can one legitimately claim to have 'weighed' the uni- 
verse or a galaxy, and even then only indirectly. Unless 
and until there is such a proof, pretending that the uni- 
verse or a galaxy has been weighed spreads ignorance, 
not knowledge. The most elementary unit of knowledge 
is a distinction made — a fact built in to every digital 
computer. Conversely, the most elementary unit of igno- 
rance is a distinction not made. 



B. Choice of a variational principle for gravity 

The step up in Sec. II from the nonrelativistic varia- 
tional princple (5j(|V(/)p + Sttkim/)) (Px = that gener- 
ates the Poisson equation = Attkii to the relativistic 
principle 5j{R— ^^n) Igl^d'h: = is an application of 
Occam's Razor, and as such can only be justified retro- 
spectively, by its consequences. Alone, this principle has, 
through its Eulcr-Lagrangc equations, little to say about 
cosmology, but when modified by inclusion of the cosmo- 
logical constant A it says something significant, namely, 
that there is gravitationally repulsive matter in the uni- 
verse and there's more of it than there is of attractive 
matter. This is a consequence that has consequences, but 
to get them all requires the next step up, to the principle 
of Eq. (6) which has in it in addition to the positive and 
negative mass densities fi and p, the scalar fields (p and ^p. 

One will recall that in deriving the field equations (7) 
and (8) only the metric components and not the densities 
or the scalar fields were varied in Eq. (6). Not varying 
the densities is normal, but not varying the scalar fields 
goes against orthodox practice. The rationale for leaving 



the scalar fields unvaried is this: In a space-time manifold 
everything is determined by the metric, and that should 
be all there is to be fixed. To vary the scalar fields would 
be to treat them as something extraneous to the metric, 
whereas their proper role should be to help describe the 
metric, not to scirvc as a source of some sort. 

Looked at another way, including the terms — -^^/i, 
— ^^/i, 2(f>~'4>.^,, and — 2 -0^ in the action integrand 
can be seen as accounting for internal constraints on the 
metric, analogously to the way including a Lagrange mul- 
tiplier A (A in the ACDM case) accounts for an external 
constraint. Just as one does not vary A (or A), one should 
not vary /x, /x, (j), or ip. In this role jj,, ft, (j), and ^ simply 
relax the field equations to allow a larger class of metrics 
to satisfy them. The relaxed field equations are as useful 
as the metrics that satisfy them, no more, no less. 

The inclusion of the scalar field ip coupled to the metric 
with polarity opposite to that of (/)'s coupling is dictated 
in the first instance by the observation that if the term 
were absent from Eq. (21), then, once R{t) surpassed 
^B/Ap (which is Rh„,,,^ when fc = 0), Eq. (21) would 
be false. Moreover, if the shoe were on the other foot 
and the term were absent, then Eq. (21) would be 
false when R{t) was less than yBjAp. Thus both (p 
and ip are needed, coupled to the geometry with opposite 
polarities. That the combination 2 4>''^(j).^ — 2 ip'^^,^ in 
the variational principle of Eq. (6) is the real part of 
2 where X •= + ^V"; suggests that space-time as 

seen here might be a restriction of a more general, at 
least partially complexified space-time geometry. 

In the drainhole model described in [10] and here in 
Sec. V only the scalar field (p was used, and it was as- 
sumed to satisfy the field equation □</> = 0, obtained 
in [10] by varying (p in the variational principle. That 
varying, however, ought to have been omitted, as = 
was already a consequence of the equations obtained from 
varying only the metric. In the drainhole model, there- 
fore, is not a source of anything, rather is just a de- 
scriptor of the geometry of space-time — and primarily, 
as seen by setting rn = in Eqs. (62) and (63), of space 
alone. That the polarity of its coupling to geometry is 
opposite to the orthodox is precisely what allows the neg- 
ative curvatures necessary to the existence of the open 
drainhole throat. To say that a scalar field coupled to 
the space-time geometry with this polarity is a 'ghost' 
or 'phantom' field associated with 'exotic matter' is sim- 
ply to indulge in misleading jargon, an error whose root 
can be traced to an uncritical acceptance of Einstein's 
assumption that energy produces gravity. 



C. Inflation and the 'big bounce' 

If the field equations dictate that R have a positive 
minimum value iJmin, then the history of the universe 
will include a bounce off that minimum; if not, then in- 
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stead of a bounce there will be a 'bang', where R = 0. 
As is evident from the definition of Pi{R) in Eq. (24), as 
well as from examining Fig. 1, for i?min to be positive it 
is necessary and sufficient that the integration constant 
B be positive (otherwise Pi{R) = —B > at i? = 0, 
which would make i?min = 0). Also evident is that the 
smaller a positive B is, the cl oser i^m in is to 0; in fact, 
as noted in Sec. Ill, iimin ~ \/B/2Ap as B 0. More- 
over, it is clear that the term 2B/3R^ in Eq. (23) is the 
term that produces rapid inflation when R is small by 
making R large, but only when B > 0. And because 
B-ApR^ - B-Ap{B/2Ap) = B/2 as B ^ 0, positivity 
of B demands that the term be present in Eq. (21), 
thus that the scalar field (j) be present in the field equa- 
tions, coupled to the space-time geometry with the un- 
conventional polarity. This is the way in which relaxing 
the field equations by including (j) allows the geometry 
to produce both a 'bounce' and the inflation that turns 
it into a 'big bounce'. Let us note that this is accom- 
plished without invoking a 'slow roll down a potential', a 
common element in many proposed inflationary scenar- 
ios. Indeed, that </> is treated merely as an aid in describ- 
ing the space-time geometry rather than as a 'physical 
source' of the geometry makes the notion of introducing 
a potential function for ^ irrelevant. 



D. The PCB model vis-a-vis the ACDM model 

After the matching — (y^p,o + A*p,o) — A^m — Mp,o in 
Sec. V the only functional difference between the expres- 
sions for in the PCB model and in the ACDM model 
lies in their terms of highest order in 1/R: Hq Vtp,/R^ in 
Eq. (65) and H"^ ^p /R^ in Eq. (66). In fully normalized 
form these equations read 

H^=Hl(nc + n^^+np^^+n^V^ (72) 



with the result that so long as $1r < 10"'"' the values com- 
puted for will match those of Table I to the precision 
shown, larger values of Hp generating larger values of 
(and younger ages for the universe). It is therefore the 
case that the PCB model is at least as consistent with 
the SNe la data as the ACDM model is. 

The omission from the PCB-model field equations of a 
radiation term proportional to 1 /i?"* goes against conven- 
tional wisdom but is not without motivation. Inclusion 
of such a term would seem to violate the basic work- 
ing hypothesis that energy is not a generator of gravity. 
Radiant energy, however, although convertible to kinetic 
energy through its influence on the motion of electri- 
cally charged objects, is not kinetic energy and therefore 
might be a gravity generator. This would presumably 
come about through a coupling of the electromagnetic 
field components to the space-time geometry. The 
difficulty is that the usual method of coupling by insert- 
ing a term proportional to F^gF^^ into the integrand of 
the variational principle has no geometrical justification 
and produces field equations of the sort that Einstein 
likened to a building with one wing made of fine marble 
(Rap — \R(lap): the othcr, of low grade wood (e. g., the 
energy-momentum tensor for F) [30] . What is needed is 
a version of the PCB model developed in the context of 
a geometrical unified field theory of the sort that Ein- 
stein spent his later years in an unsuccessful quest for. 
Although I have been able to construct such a theory as 
a hybrid of Weyl's conformal theory and Kaluza's five- 
dimensional theory [31, 32], the manner in which that 
theory's variational principle and field equations connect 
the electromagnetic field to geometry does not suggest a 
way to relate F to the space-time geometry in the present 
context. For this reason I have chosen to omit the radi- 
ation term from the field equations governing the PCB 
model. Such a term can, of course, be reinstated on an 
ad hoc basis if a sufficiently compelling reason for doing 
so presents itself. 



(which is Eq. (42)) and 

H^ = Hl[n^ + n^^+nM^+HK-^}j- (73) 

Absent their highest-order terms, and with the obvious 
identifications among the remaining ils, these equations 
would produce, by way of the corresponding formulas for 
the distance modulus /x, exactly the same best fits to the 
SNe la data, the parameters and the statistic being 
those of Table I. The rationale behind the construction 
of Table II was that, because the SNe la data extend 
into the past only to 2: = 1.755, restoring the fie term in 
Eq. (72) while retaining the Ob = best-fit values of the 
other parameters will yield a fit with negibly different x^- 
The same reasoning could be expected to apply also to 
restoring the Or term in Eq. (73). This I have tested, 



E. Dark matter, dark 'energy', and the 
'Cosmological Constant Problem' 

In the ACDM model dark matter is lumped with bary- 
onic matter in Om as primordial matter whose density is 
proportional to 1/R^. The unseen, unknown stuff in Oa, 
whose density has no R dependence, is called dark 'en- 
ergy' in deference to Einstein and his assumption that 
energy is a source of gravity. In the PCB model only 
baryonic matter is treated as belonging to the primordial 
Op sector with 1 /R^ density dependence, dark matter be- 
ing shifted to the Oq sector with density kept constant 
by continuous creation of the dark C-matter. And there, 
instead of dark 'energy' represented by a cosmological 
constant A many orders of magnitude too small to be 
consistent with various origins proposed for this energy. 
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one finds the gravitationally attractive dark matter parti- 
cles accompanied by their gravitationally repulsive back 
sides, their net effect represented by the constant Ac of 
Eq. (20), functionally equivalent to A and of the requi- 
site smallness. This interpretation of 'dark energy' as the 
back side of dark matter provides a self-consistent solu- 
tion of the so-called 'Cosmological Constant Problem'. 



F. Ratio of dark matter to baryonic matter 

Because in the ACDM model the densities of dark mat- 
ter and baryonic matter are both proportional to l/R^, 
the ratio of their densities is constant in time. Estimates 
of that ratio obtained from observations of visible struc- 
tures such as galaxies and galactic clusters and super- 
clusters as they exist at the present epoch in the vicinity 
of our galaxy are therefore used in studies of the forma- 
tion of these structures in the distant past. In the PCB 
model, on the other hand, the dark matter density is held 
constant through continuous creation of C-particles, so 
the ratio /ic I MP of dark C-matter density to baryonic P- 
matter density is proportional to . Consequently, the 
value of this ratio at the epoch corresponding to redshift 
z is 7(z) = 7o/(l + 2:)^, where 70 is the ratio at present, 
several possible values of which are inchided in Table III 
for comparison purposes. Taking this dependence of the 
ratio on redshift as real would no doubt influence the 
outcomes of studies of structure formation. 



G. Continuous creation of C-matter tunnels 

The stipulation that the net density /xc -l- of 
C-matter stay fixed while the universe is expanding re- 
quires that the drainhole tunnels considered to be the 
particles of C-matter come continuously into existence. 
By what mechanism might this happen? Wheeler's no- 
tion of a 'quantum foam' of wormholes popping into and 
out of existence would not provide a satisfactory explana- 
tion, for to keep the density constant each newly created 
tunnel must either remain in existence or else upon dying 
be multiply replaced by new ones. A more useful idea is 
that the tunnels arise from the stretching of space as the 
universe expands. Once started, this would be a regen- 
erative process, each new tunnel causing more stretching 
that would generate additional tunnels. Such a process, 
moreover, could be expected to produce particles in a 
given region of space at a rate proportional to the rate of 
increase of the volume of that region, thus maintaining a 
constant particle density. If, as was assumed in Sec. V, 
these particles are all created with the same mass and 
size parameters mc and nc, the net density would be 
held constant. As to how the process could begin, one 
looks to the primordial P-matter tunnels, presumed to 
have always been present, pushing the universe toward 



expansion. In the postbounce era the C-particles would 
add their contribution, the ultimate result being an ac- 
celeration of the expansion. In the prebounce era the C- 
particles, abundant in the distant past, would be dying 
out with the shrinking of space as the universe contracted 
into the bounce. 

If the primordial P-matter has 'always' existed, then 
the question of what 'caused' it to exist is better left 
to the philosophers. Conceivably, however, the P-matter 
tunnels came into existence as dark matter particles the 
same way the C-matter tunnels did, but in the inflation- 
ary space-stretching era after the bounce. Only later 
would their entrance portals have acquired the decora- 
tions that turned them into baryonic matter. The re- 
ally fundamental unanswered question, though, is how 
did space and time come to exist, for, time excepted, 
all the elements of the PCB model — the tunnels and 
the scalar flelds — are merely features of the geome- 
try of space. In fact, the question can be reduced to 
the existence of space alone, for the Kaluza-Weyl the- 
ory described in [31, 32] is based on a construction that 
produces time (thus space-time) from three-dimensional 
space and, repeated, produces a secondary time (thus 
space-time-time) from space-time. And in the space-time 
stage of that theory there occurs naturally a drainhole 
model in an expanding universe which repels on the low 
side more strongly that it attracts on the high side and 
needs no auxiliary scalar fleld for its description [33]. 

That all the C-particles are created with the same mass 
and size parameters is only a convenient assumption. If 
this assumption is relaxed, then the densities /xc and /zc 
might vary, in which case the proportionality of \xc + 
/xc to 1/i?^ could cease to be exact. This variability 
could then be used in an ai hoc way to produce some of 
the putative behaviors that proposed modifications of the 
ACDM model put forth under the rubric 'quintessence' 
are intended to allow. 



H. Voids, walls, filaments, and nodes 

The PCB model, just as the ACDM model, treats 
matter as uniformly distributed in space, an apparently 
reasonable assumption on a large enough scale. Obser- 
vationally, however, a major fraction of space appears 
to comprise voids nearly empty of matter, bordered by 
walls, filaments, and nodes in which resides most of the 
matter, both the dark and the visible (baryonic). The 
formation of this cellular structure is generally consid- 
ered to be an evolutionary product of early small fluctu- 
ations in the density of primordial matter. Explanations 
along that line are attempts to realize what Peebles has 
called his "perhaps desperate idea . . . that the voids have 
been emptied by the growth of holes in the mass distribu- 
tion" [34] . If besides matter that attracts gravitationally 
there is also repulsive matter, then that idea can be rein- 
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terpreted and the desperation perhaps reversed, by the 
simple observation that, whereas attractive matter wants 
to congregate, repulsive matter is reclusive, pushing all 
matter, attractive or repulsive, away from itself. 

In the PCB model each tunnel particle would presum- 
ably be created with its entrance and its exit close to one 
another in the ambient space (if far apart, all the better). 
In the ambient space the entrance would attract the exit, 
but the exit would repel the entrance more strongly, so 
the two portals would drift apart. Apply this to a mul- 
titude of such particles and you would expect to see the 
exits spread themselves over regions from which they had 
expelled the entrances, regions therefore devoid of ordi- 
nary (i. e., attractive) matter. The entrances, on the 
other hand, being brought together by both their mu- 
tual attractions and the repulsion from the exits, would 
aggregate into walls, filaments, and nodes on the bound- 
aries of the void regions, just as observed in the real uni- 
verse. What is more, the walls, filaments, and nodes so 
created would likely be, in agreement with observation, 
more compacted than they would have been if formed by 
gravitational attraction alone, for the repulsive matter in 
the voids would increase the compaction by pushing in 
on the clumps of attractive matter from many directions 
with a nonkinetic, positive pressure produced by repulsive 
gravity, a pressure not to be confused with the negative 
pseudo-pressure conjectured in the confines of Einstein's 
assumption to be a producer of repulsive gravity. In this 
way the PCB model explains qualitatively the observed 
cellular structure of the universe with minimal or no re- 
liance on fluctuations in the primordial matter density. 

The validity of this explanation of the cellular struc- 
ture could be tested independently of the tunnel parti- 
cle hypothesis by numerical simulation of evolution of 
an expanding universe with both attractive matter and 
repulsive matter distributed throughout. An interesting 
topic for investigation would be the propagation of light 
through the voids, particularly as the repulsive matter in 
a void might tend to distribute itself nonuniformly. 



I. Geometric units and inertial mass 

In Sec. V, in order to populate Table III, the radius 
r(2mc) of the C-particle throats was set equal to vari- 
ous multiples of the proton- neutron inertial mass mp.n, 
converted to geometric units to give the radius its req- 
uisite units of length. Because the conversion from in- 
ertial mass units to geometric units (Ig — )• G/c^ x Ig « 
7.426 X 10~^^ cm) involves Newton's gravitational con- 
stant G, such a conversion would seem to be at odds 
with the assumption that inertial mass does not gener- 
ate gravity. The resolution of this apparent inconsistency, 
dictated by the proposition that everything that exists is 
one or another aspect of the geometry of the universe, is 
that the gram is an artificial unit of inertial mass used 
for convenience, and that the proper unit of inertial mass 



is a length, 7.426 x 10~^^ cm, for example. That neces- 
sitates that G/c^ = 1. The 'geometry is all' proposition 
dictates also that the second must be considered an arti- 
ficial unit standing in for a length unit. Taking that unit 
to be the light-second makes c = 1 and, consequently, 
G = 1, which equations are commonly taken as the defin- 
ing equations of geometric units. Considering the proper 
units of inertial mass and time to be units of length does 
not, of course, prevent the use of the gram and the second 
in their customary contexts. 



J. Galactic nuclei as drainholes 

If instead of a blackhole at the center of our galaxy 
there is a drainhole, what might it look like? Prom 
afar the behavior of matter and radiation around such 
a drainhole would not differ in appearance from the be- 
havior of matter and radiation around a blackhole of the 
same mass and approximately the same size. For both 
a Schwarzschild blackhole of mass m and a drainhole of 
mass m the radial equation of motion of a test particle is 



dr'^ 



m 



,(^^+(p-3m)(^-J 



(74) 



with r(/7) given by Eq. (63) for the drainhole and r(p) = p 
for the blackhole. In both cases r is the proper time 
of the test particle. Equation (74) shows that for the 
drainhole as well as for the blackhole 2>m is its 'star- 
grabbing' radius, as any star whose orbit takes it down to 
p = Srn or below will never again have dp/dr > 0. The 
shared metric of the cross sections of constant t being 
dp^ + [p) dil^ , the cross sections of the blackhole are 
euclidean, so the geodesic radius 3m of the sphere from 
which no star can escape is also its areal radius. For the 
drainhole the situation is different: the areal radius is 
r(3m), which by increasing the size parameter n can be 
made as large as you like. 

For purely radial motion the acceleration —m./r'^{p) 
in Eq. (74) is the (nonvectorial) gradient [p /2)'{p) of 
the specific kinetic energy ^tP' of observers free- falling 
from rest at p = oo. It attains its maximum strength 
m/r^(2m) at the choke point of the drainhole throat. A 
measure of the tidal effects in the radial direction is the 
tidal gradient, given by 

.<n\. (^^^(\ 2(p-2m) m 

\2 J [p — my + n-^ — m'^ r'^KP) 

which vanishes at the choke point and reaches a positive, 
stretching maximum 



^/37 



2n + V^m r'^{p+) 



at p+ = 2m H — -= 



n 



(76) 
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and a negative, compressing minimum 

-\/3m 1 r, 

(77) 

For a fixed m, Smax is a monotonia function of n, decreas- 
ing from 

3V3(2-V3)2e-(3-v^) 
lim s„,ax = 78 

to 0, while Smin increases monotonically from 



lim 



^min — 



3(12 + 7x/3)e-(3+V3) 
2m? 



(79) 



to 0, as n increases from m to oo. 

Observations of stars near the galactic center put the 
central supermass constraining their motions at mgc ~ 
4.31 X 10^ M© = 6.36 x 10^ km [35]. In ordinary units 
1/mgc = cV(6.36 X 10^ km) = 1.41 x lO'' km/s^ and 
l/ml^ = cV(6.36 X 10^ km)^ = 2.22 x lO^^ (km/s2)/km. 
Building on these, with the ratio njm^c as an input 
parameter, Tabic IV exhibits the corresponding values 
of the drainhole areal radius r(2m), the star capture 
areal radius r(3m), the maximum radial acceleration 
m/r^(2m), and the tidal gradients Smax and Smin- 

Stars, gases, and light falling into a supermassive 
blackhole encounter unbounded tidal stretching that is 
said to strip them of their identities as they approach 
and reach the central singularity, never to be seen again. 
Falling into a supermassive drainhole gases and light (and 
even stars if the ratio n/mgc is large) can retain their 
identities and exit the drainhole on the other, gravita- 
tionally repelling side. This brings the following ques- 
tion: If supermassive drainholes at galactic centers are 
commonplace, and their exits reside in our (part of the) 
universe, where should we expect to see those exits, and 
how might we identify them? In accordance with the no- 
tion that the exits of drainhole tunnels tend to isolate 
themselves by pushing all other matter away, the obvi- 
ous places to look would be the great voids filled with the 
exits of dark C-matter tunnel. Presumably they would 
appear as compact objects surrounded by a halo of out- 
flowing gas, their light intrinsically blueshifted from its 
journey downward in the drainholes' gravitational poten- 
tials. The brightnesses of such objects would depend on 
the rate at which the gases and light were sucked into the 
drainhole entrances. If very bright, they might appear as 
isolated quasars. If dim, they might be present in the 
voids but unseen. 

With a little imagination one can see how a supermas- 
sive drainhole might arise from a plenitude of baryonic 
P-matter and dark C-matter drainhole tunnels with their 
entrances brought tightly together by their mutual at- 
tractions. Begin with a pair of them, T\ and Ti say. If 



not only their entrances but also their exits were close 
together, then the coalescing of T\ and Ti into a single 
tunnel might be feasible. If, however, as I have proposed, 
their exits had receded far from one another, then the 
only way T\ and Ti might reasonably coalesce into one 
tunnel would be to join at their entrances and close off 
into a single tunnel connecting their original exits. If 
that occurred, then the 'gravitational ether' flowing into 
their entrances would be diverted to the nearby entrances 
of one or more tunnels, thereby increasing their masses, 
both the attractive and the repulsive. The remnant of 
T\ and T2, no longer connected to the ether flow, would 
soon become a sterile, gravityless tunnel connecting two 
distant places out in the voids. After many repetitions 
of this process there would be left a single tunnel large 
enough to accommodate all the combined ether flow of 
the original P- and C-matter tunnels — a supermassive 
drainhole. 

One can equally well imagine that all the tunnels as- 
sociated with baryonic matter and dark matter concen- 
trated at a galactic center would retain their identities 
and their gravity, both the attractive and the repulsive, 
no matter how tightly packed or rapidly moving their 
entrances might be. At the tunnel exits there would be 
nothing new to look for, all the attractive gravity concen- 
trated at their entrances being dispersed in the voids as 
negative gravity in the same way as it was before the con- 
centration took place. Here there would be no merging 
of tunnels, thus no altering of the topology of space. Left 
to be analyzed mathematically would be only the simpler 
changes of spatial topology associated with the creation 
and destruction of the individual drainhole tunnels. 



VII. SUMMARY 

Cosmic inflation, deceleration, acceleration, dark mat- 
ter, dark 'energy', and a 'big bounce', not to mention 
coasting, voids, walls, fllaments, and nodes, have been 
found wrapped in one coherent package, namely, the 
variational principle (5/[i2 — ^^{/J. + jl) + 2 0''''(/).-y — 

2^-'''V.7] bl^c?^ = and the field equations it implies. 

Arriving at this variational principle depended on giving 
up Einstein's assumption that inertial mass and energy 
produce gravity. 



VIII. NOTE 

This is a revision and major extension of the earlier 
version(s) [36], retitled to reflect the revision. 
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TABLE IV: For a drainhole of mass parameter m = mgc and size parameter n, where mgc ~ 4.31 x 10 Mq = 6.36 x 10 km 

in geometric units (the current best estimate of the supcrmass at the center of our galaxy), and a range of values of the ratio 
n/m-gc, this table exhibits the corresponding values of the drainhole areal radius r(2m), the star capture areal radius r(3m), 
the maximum radial acceleration m/r^{2m), and the maximum and minimum tidal gradients Smax and Smin- 



Inputs 


Outputs 


Tt 

mgc 


r(2m) (km) 


r(3m) (km) 


m n -2^ 

-r-;- r (km S 

r2(2m) ^ ' 


Smax (km 


s-2 


/km) 


Smin (km s ^/km) 


1 


1.73 X lO'^ 


2.10 X 10^ 


1.91 X 10^ 


1.17 X 


10" 


-4 


-7.08 X 10" 


-4 


10 


7.37 X lO'' 


7.41 X lO'' 


1.05 X 10^ 


1.01 X 


10" 


-6 


-1.15 X 10" 


-6 


100 


6.46 X 10* 


6.46 X 10* 


1.37 


1.39 X 


10" 


-9 


-1.41 X 10" 


-9 


1000 


6.37 X 10^ 


6.37 X 10^ 


1.41 X 10"^ 


1.44 X 


10" 


12 


-1.44 X 10" 


12 


10000 


6.36 X 10^° 


6.36 X 10^° 


1.41 X 10"* 


1.44 X 


10" 


16 


-1.44 X 10" 


15 


100000 


6.36 X 10" 


6.36 X 10" 


1.41 X 10"® 


1.44 X 


10" 


18 


-1.44 X 10" 
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